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INTRODUCTION 


DESPITE THE LARGE NUMBER of studies on the feeding 
biology of deposit-feeding bivalves, questions concerning 
the dependence of growth upon sediment composition 
and the application of allometric relations for describing 
growth remain to be answered. To better understand the 
nutritional basis for growth of deposit-feeders, HYLLEBERG 
& Gattucct (1975) studied the feeding selectivity and 
growth of Macoma nasuta (Conrad, 1837), a deposit- 
feeding intertidal bivalve of the west coast of North 
America (Pouro, 1973). 

NEWELL (1965) worked with Macoma balthica (Lin- 
naeus, 1758) and demonstrated that the animals’ ener- 
getics are highly dependent upon the microorganisms 
which coat the surface of sand particles. Macoma nasuta 
is also capable of removing microorganisms from particle 
surfaces. A silt-clay environment is considered to be nutri- 
tionally richer because a given weight of fine silt-clay 
would have a greater surface area and thus microorganism 
content than that of a coarse sand. In addition to micro- 
organisms, detritus is selected for and makes, at its early 
stage of decomposition, some contribution to the nutrition 
of M. nasuta (HyLLEBERG & Gattucci, 1975). The sig- 
nificance of suspension feeding is incompletely known in 
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tellinaceans, although the possibility is often mentioned, 
e. g., in BRAFIELD & Newer (1961). In M. nasuta the 
significance of suspension feeding was not determined 
(HyLLEBERG & GALLUCCI, op. cit.). 

There is an extensive literature on growth in bivalves, 
some of which is quantified in the form of various types 
of mathematical growth models. Analytical growth 
models usually entail at least 2 parameters, a maxi- 
mum body size and a growth constant. If growth is 
noted by a linear measurement (e. g., length) the maxi- 
mum size would be Le and the corresponding weight Wa. 
Growth rates in time could then be di(t)/dt and dw(t)/dt 
and the growth constant would be k. The well known von 
Bertalanffy model uses maximum sizes, growth rate and 
a third parameter te (von BERTALANFFY, 1938) and other 
more elaborate models contain additional parameters 
(SourHwarp & CHapMaAN, 1965; FLETCHER, 1974). 
However, the estimation and testing of parameters must 
be done with care to avoid problems associated with the 
fact that parameters may not be independent of each 
other (GaLLuccI & QuINN, in manuscript). 

This paper reports on comparative growth studies of 
Macoma nasuta in silt-clay and medium coarse sand at 
approximately the 0.0 tide line. Concurrent studies were 
done in the laboratory (HyLLEBERG & GaLLucct, 1975) 
where the proportion of the contribution to growth during 
different modes of feeding could be controlled. In contrast 
to analytical growth models, time (t) is not the independ- 
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ent parameter. Instead, the data are reported as weight 
(ash free) w(t) as a function of length /(t). The empir- 
ical, allometric model of the form w = al” shows that the 
maximum length is a critical parameter, as is the growth 
rate q = dw/dl. 


METHODS 


The field experiments were done in Garrison Bay, a 
small bay of San Juan Island in the San Juan Archipelago, 
Washington (Figure 1). The bay is approximately 3km 


open area 
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Figure 1 


Garrison Bay, San Juan Island, Washington, USA 


long and 0.8km wide with a classic sediment profile: fine 
silt and clay substrate at the closed end where the anaero- 
bic layer is very near the surface and coarse sand at the 
open end. The experiment to evaluate growth consisted of 
burying 9 measured Macoma-nasuta in the open area of 
the bay and 9 in the closed area. In each set of 9, 3 were 
juvenile (pre-reproductive), 3 were in the mid-size range 


and 3 had sizes near Lx. The 9 clams in the closed area 
were selected from about 60 that were dug until the 
appropriate sizes were found. They were all dug from a 
region adjacent to where the 9 were buried. A similar ex- 
planation applies to the open area. Each clam’s length 
was measured and the clam was placed in a plastic 1-liter 
capacity container which was cut at the sides and bottom 
and covered with a 1mm window screen to allow water 


Figure 2 


Data Points and the Curve w = a |b from the Open Area 
Field Experiment (FO) 
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Figure 3 


Data Points and the Curve w = a lt from the Closed Area 
Field Experiment (FC) 


transport across the boundary. Each container was filled 
with sediment from an adjacent location and buried flush 
to the sediment. All of the 9 containers were buried ad- 
jacent to each other and parallel to the water line in 
areas that visually appeared similar with respect to living 
and decaying algal cover (Enteromorpha intestinalis and 


Ulvoids). Therefore, differences in exposure to organics 
would largely be a consequence of sediment differences. 
Since sediment in the containers was sieved to remove 
possible competitors, any competition was due to epifaunal 
colonizers such as Transennella tantilla (Gould, 1853). 
After 40 days of growth, starting in mid-July, 1974, the 
clams were dug up and the lengths measured. Shell length, 
the longest possible postero-anterior distance, was meas- 
ured to within 0.1 mm with a vernier caliper. The ashfree 
weights of the animals that had been growing for 40 days 
were found by drying animal tissue for at least 24 hours 
at 80°C and burning at 550°C for 24 hours. It was de- 
termined that M. nasuta valves contain no more than 5 
mg of organic matter. Figures 2 and 3 contain the data 
points from field experiments on weight w(t) versus 
length /(t) axes for the open and closed areas, respectively. 

The differences in the particle size composition of the 
growing regions were quantified. Two sediment cores each 
were taken from the open and closed areas of the bay. 
The cores were dried at 100°C for 4 hours; the dried 
sediment was picked free of large organics, put on an 
automatic shaker and passed through a series of sieves of 
mesh sizes: 1.980, 0.495, 0.246, 0.124, 0.063mm. The 
sediment retained by each sieve and that passed through 
the 0.063 was weighed to get the percent of sediment in 
each size range. The averaged cumulative percentage 
weight curves of the 2 sediments are similar in shape but 
the closed end area always contains a silt-clay component 
at least 2 times larger by weight than the open end. 

The differences in the organic content of the 2 areas 
was found by loss on burning. Organic content of the 
sediment was determined by taking a sediment core from 
each area of the bay, freezing it for 4 days, removing the 
top 4cm from the cores, drying and burning as above. Loss 
on burning was calculated as loss/mg dry weight of sedi- 
ment. Sediment from the closed area contained twice as 
much combustible organics assediment from the open area. 
Presumably much of the difference is due to there being 
a greater number of microorganisms and detritus in the 
sediment from the closed area of the bay. 

Salinity in the bay is relatively constant at about 28 to 
30%. during the primary growing time, but at specific 
“run-off areas” the salinity of the water that fills a hole 
dug in the intertidal zone may be as low as 16%,. Neither 
growth study was in such an area. Water temperature in 
the open area is about 12°C and on warm sunny days it 
is about 2°C higher in the closed area. 

In late August an estimate of the average bivalve 
density (Macoma nasuta, M. inquinata, M. inconspicua 
(balthica), Protothaca staminea, and other large suspen- 
sion feeders) in each area of the bay in the vicinity of the 
growing site was found by digging a 1.5m x 25cm x 60cm 
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trench and sieving all sediment through a mesh size of 6 
mm, thus recovering most of the small size classes. The 
M. nasuta density in the open and closed areas was ap- 
proximately the same and, although the size-frequency 
distributions (Figure 4) are translated horizontally, they 
are similar. 
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Figure 4 


a. Size-frequency Curve of Macoma nasuta from the Open Area 
b. Size-frequency Curve of Macoma nasuta from the Closed Area 


Predicted growth curves of the form 

w =a? (1) 
were fitted to the field data plotted in Figures 2 and 3 and 
are the solid curves drawn in these figures. The allometric 
parameters, a and b, of the curves are in Table 1. These 
parameters for the field experiments in the open (FO) 
and the closed areas (FC) were statistically compared to 
each other. Similar statistical- comparisons of parameters 
a and b were done with the growth data from compatible 
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Table 1 


The values of allometric parameters and the 
corresponding asymptotic variances from studies 


done in the field in the open (FO) and closed (FC) areas 


Parameter Asymptotic Parameter Asymptotic 
Area a variance variance 
Open (FO) 0.014084 0.010117 2.8655 0.19233 
Closed (FC) 0.30845 0.32049 2.0495 0.27526 


laboratory studies (HyLLEBERG & Garucci, 1975) and 
are given in Table 2. The laboratory experiments were 
conducted simultaneously. The data from these studies 
are plotted on weight versus length axes in Figure 5 
(LO) and Figure 6 (LC) which correspond to labora- 
tory simulations of growth in the open and closed areas, 
respectively. The solid curves correspond to equation (1) 
with the estimated values of a and b. 


Table 2 


The values of allometric parameters and the 
corresponding asymptotic variances from studies done 
in the laboratory to simulate growth in the open (LO) 

and closed (LC) areas 


Parameter Asymptotic Parameter Asymptotic 
Area a variance variance 
Open (LO) 0.06388 0.073293 2.4654 0.30307 
Closed (LC) 1.4604 4.4586 1.5948 0.82595 


The allometric parameters a and b for the laboratory 
and field experiments were also compared. The test of 
comparison is Rao’s Chi-Square Test (Rao, 1973: 389 
© S 

The non-linear equation (1) was fit using the least 
squares BMD X85 program on a CDC 6400. This par- 
ticular equation was used for 2 reasons: the experimental 
data are clearly of this form with no apparent inflection 
or turning points and because the equation, known as the 
allometric relation, has found wide application in growth 
studies of many organisms (BEVERTON & Hott, 1957; 
ALDRICH, 1974). 
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Figure 5 


Data Points and the Curve w = a lb from the Open Area 
Laboratory Experiment (LO) 


It is well known that the parameters a and b may be 
estimated from a log - log transformation of (1) and lin- 
ear regression. Parameter a would correspond to the verti- 


Figure 6 


Data Points and the Curve w = a lt from the Closed Area 
Laboratory Experiment (LC) 


cal intercept and b to the slope. This method of estimating 
a and b is common in early work but is now unacceptable 
due to the bias which may be introduced by the trans- 
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formation (Grass, 1967; Zar, 1968). Non-linear regres- 
sion is preferable in this age of calculators and computers. 
Nevertheless, the transformation is a useful way to demon- 
strate the form of, and to compress, the data. 


RESULTS 


The growth data from the field experiments are in Fig- 
ures 2 and 3. If the solid curves calculated from (1) are 
superimposed an intersection is observed. The point of 
intersection is the value / = l, which satisfies (see Table 1) 
0.308)" == 0.014/7°* 

which is 1; = 43.94mm. That is, for values of l(t) 
which are less than J, the biomass w(/) at any allowable 
l-value is greater in the closed area, where the sediment 
is organically rich silt-clay, than in the open area. How- 
ever, the presence of an intersection implies that the rate 
of growth of biomass per unit length, dw/dl, must be 
greater in the open area for l < l. If L»° is the maximum 
size of open area clams then, for l < l< L°, the w(l) 
of open area clams is greater than that of closed area 
clams. If L«* is the maximum size of closed area clams it 
is apparent from Figures 2, 3, 4 that Le° > Le«®. Thus 
dw/di is not simply a function of sediment type (or equi- 
valently, food availability) but must be a function of food 
availability and shell size, also. 

The derivate of w(t) with respect to l(t) isq. From (1) 

q — dw/dl Ja bi 
Thus in the open area q=0.0401' ” and 
in the closed area q = 0.6311% (see Table 1). 

Plots of q versus / for each area show that the curves 
intersect and that the rate of change of biomass with 
respect to length is almost linear in the closed area and 
almost quadratic in the open area. That is, the instant- 
aneous rate of change functions have the same form we 
attribute to them from arguments based on finite rates of 
changes in Figures 2 and 3. There is, of course, no need 
for the curves of w versus l to cross at the same value as 
the curves of q versus l. 

Curves fit from laboratory experiments (HYLLEBERG & 
Garucci, 1975) using the same sediment types show the 
same qualitative behavior. The intersection of the curves 
in Figures 5 and 6 occurs at the | = l which satisfies (see 
Table 2) 

0.064174 — 1.460139 
which is /2 == 36.40mm. Thus the laboratory studies show 
that closed area clams have a greater biomass than open 
area clams for l < ls. This result is compatible with, but 
less than the values found in the field study. 

The comparison of the growth in different sediments 
in the field and in the laboratory is of more value when a 
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quantitative statement can be made about the probabili- 
ty of these differences in parameter estimates occurring 
again. In Table 3 the results of a Chi-square test (Rao, 
1973: 389 - 391) are given and the estimates compared 
at the 99 and 95 percent probability levels. 


Table 3 


FO, FC, LO, LC are field study open area. field study 
closed area, lab study open area, lab study closed area. 
Rao Chi-square test is applied to estimates of parameters 
a and b using asymptotic variances. Parameters 
significantly different in the locations indicated at 
a = 0.05 are indicated by * and those different at 
a = 0.01 are indicated by **. 

Chi square values x?o.0o5 = 3.84, x?o.01 = 6.64 


Null hypothesis Parameter Parameter 
HS a b 
FO = FC 5.90* 44.59"* 
LO = LC 0.88 8.81** 
FO =LO 4.08* 11.18** 
IC = ILC 0.60 2.39 


The size frequency distributions of Macoma nasuta 
(Figure 4) taken in late August from each of the 2 areas 
show a tendency for larger (longer) animals to occur 
more frequently in the closed area where the larger peak 
is shifted to the right by about 3mm. The total number 
of these deposit feeding animals is, however, approxi- 
mately the same in both areas. 


Table 4 


Average change in length A/ over time 
for clams in field experiments 
in the open (FO) and closed (FC) areas 


Size range 
Area smallest mid largest 
Closed (FC) 4.1 0.40 0.15 
Open (FO) 2.7 0.70 0.60 
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The data on changes in length (Al) over the 40 day 
growing period in the field are given in Table 4 and show 
unambiguously that there is a greater growth in the 
closed area in the smallest size classes. Furthermore, the 
data from the largest size class, all clams which were 
greater than 43mm, increased more over time in the open 
area than those in the closed area. 

We have used L= to discuss the maximum size attain- 
able in a sediment type. However, we have not calculated 
values from the data because, by definition, L» corre- 
sponds to the value (in the von Bertalanffy model) where 
di/dt = 0, and reliable /(t) values are difficult to de- 
termine. 


DISCUSSION 


The data from both the field and laboratory studies 
are noisy in the mid-size ranges and less so at the extreme 
sizes. Furthermore, data from the richer closed, silt-clay, 
area generally have the greater variability which is de- 
monstrated in Tables 1 and 2 from the larger asymptotic 
variances in the parameter estimates of growth in the 
closed area. It can also be shown (Hylleberg, unpub- 
lished data) that cleaned valves collected in the silt en- 
vironment have a much greater variability about a re- 
gression curve of shell weight on length than valves from 
the coarser sand area. Increased variability in growth data 
may be related to the presence of increased numbers of 
some deposit feeding species who biogenically rework the 
sediment and some suspension feeders who modify the sur- 
face sediment by their presence (RHODES & Youne, 1970; 
ScHERBA & GatLucct, 1976). The results of the latter 
Study suggest that a sheltered silt-clay environment may 
act as a buffer against the extremes of seasonal fluctua- 
tions in abiotic variables. It is probable that such areas are 
more heterogeneous and thus contain a greater number of 
patches for deposit feeding animals. A possible conse- 
quence is an increase in the variability of growth data 
with the potential for generally larger growth. The fact 
that the curves in Figures 2, 3 and 5, 6 are qualitatively 
similar with the intersections occurring in comparable 
Size ranges is further evidence that the reported growth 
characteristics are not artifacts. 

The exponent b in both laboratory and field studies is 
Significantly different at œ = 0.01 when values for the 
open and closed areas are compared. Although the co- 
efficient a in the laboratory study is not significantly 
different at either œ level (in comparison of open and 
closed areas), a significant difference at the a= 0.05 level 
is found between the parameter values estimated from the 
2 areas in the field study (Table 3). The results from tests 


of significance are certainly influenced by the above noted 
higher asymptotic variances associated with the parame- 
ter a, evaluated in the closed area in the field and labor- 
atory data. 

The data from the histograms of size-frequency (Fig- 
ure 4) indicate that more Macoma nasuta in the closed 
area are in the largest size classes. Since the number of 
M. nasuta seems to be comparable in the 2 areas but the 
closed area animals generally have a greater biomass per 
shell length, intra-species competition for food cannot be 
eliminated. Similarly, inter-species competition among 
polychaetes and deposit-feeding bivalves other than M. 
nasuta cannot be eliminated. There is no evidence of pre- 
dation on larger sized M. nasuta. There is, however, the 
possibility of predation on animalsin the 0-year class or of 
an indirect mortality such as trophic group amensalism 
(Ruopes & Younc, 1970) which would not reveal itself 
via large amounts of small, empty valves. These latter 
processes are considered as probable factors in the occur- 
rence or non-occurrence of species in a patch (GALLuccI, 
in preparation). Since our calculations indicate differen- 
ces in the abundance of food items in the 2 areas, 7. e., 
predominance of diatoms in the open, versus bacteria and 
detritus in the closed area, we are inclined to suggest this 
as an explanation for the observed growth behavior. HYL- 
LEBERG & GaxLuccr (1975) tried to confirm this sug- 
gestion by a simultaneous laboratory study where the 
growth value of food sources from the open and closed 
areas could be analyzed. However, problems associated 
with selective feeding by sorting in the mantle cavity pro- 
vided less than definitive results. In fact, MAURER (1969) 
reports a bivalve-length-sediment study in which the sizes 
of different species are positively correlated, negatively 
correlated and uncorrelated with the sediment particle 
size distribution. 

We show that deposit feeding Macoma nasuta have 
more biomass and add biomass/unit length more slowly 
in the organically richer closed area, for lengths less than 
about 43mm. This same population will grow larger on 
the average than open area clams and, for lengths greater 
than about 43 mm, biomass per unit length is added faster 
than for the open area clams. The data in Table 4 also 
show a decreased rate (with respect to time) of growth for 
large clams from the closed area as compared to the open 
area. These results indicate an environmentally related 
dependence between maximum size and dw/dl. Growth 
parameters L- and We are often estimated from each 
other using the allometric relation and Le is often estima- 
ted by fitting the von Bertalanffy growth law (BEVERTON 
& Hort, 1957). However, it does not necessarily follow 
that dw/d/ and di/dt are always related to each other in 
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a straightforward way analogous to the mathematical re- 
lationship of the chain rule derivative of w(t) with respect 
to t. Despite this, comparisons are made of populations’ 
growth or of environments’ capacity for growth using 
some of these quantifiable parameters and calculating 
the others. Our results indicate that, at least for deposit 
feeding bivalves, the biological relationships may be more 
complicated. The inter-relationships suggested here are 
usually masked in biomass versus time or length versus 
time studies, or both. 

The allometric relation (1) between biomass and length 
is a natural function to use to describe such data. The 
exponent values are between 2 and 3 as expected from 
arguments based on the relationships between biomass, 
volume and linear dimensions and from previous studies 
on marine fin fish (BEvERTON & Hott, 1957). Further- 
more, despite variability in the mid-range, the data also 
have this polynomial form. 

The estimated parameters a and b are, subject to the 
preceding, indicators of the conditions of growth the dif 
ferent populations experience. If a and b could be esti- 
mated from growth over a single year, anomalies of that 
year’s growing conditions might govern the estimates. 
However, a and b are estimated from clams growing from 
1 to about 6 years, so anomalies tend to cancel (trends 
would not cancel). Since Macoma nasuta is a dominant 
animal in the mud flat community of Garrison Bay ( GAL- 
LUCCI, in preparation) the parameters a and b may be 
considered as community indices whose stability is quanti- 
fiable. 


SUMMARY 


In this report we demonstrate the applicability of 
the allometric relation between biomass and length in 
Macoma nasuta, a deposit feeding bivalve, and statistical- 
ly test the null hypothesis of no difference in biomass ver- 
sus length curves for a population growing in a rich silt- 
clay environment and a population growing in a medium- 
coarse sand environment. The allometric parameters are 
shown to be statistically different, establishing that there 
is a difference in the way biomass increases with respect 
to valve length in each environment. 

The curves of biomass versus length in the 2 environ- 
ments intersect at 43mm. Biomass values in the richer 
environment are greater for lengths less than 43 mm. The 
biomass values are greater in the sand environment at 
lengths between 43mm and the maximum length (ap- 
proximately 50mm). Rates of increase of biomass with 


respect to length are indicated. A size-frequency study 
demonstrates that there is a greater number of larger bi- 
valves in the richer environment. 

Simultaneous studies done in the laboratory (HYLLE- 
BERG & GaLiucct, 1975) yield similar growth curves 
and generally confirm the field results. 

Both field and laboratory results point to increased vari- 
ability in silty environments which seems to translate to a 
greater heterogeneity with meaningful implications for 
understanding sediment inhabiting animals. 
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